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We study the mechanical properties of semi-flexible polymers, when the contour length of the 
polymer is comparable to its persistence length. We compute the exact average end-to-end distance 
pq: and shape of the polymer for different boundary conditions, and show that boundary effects can 

lead to significant deviations from the well-known long-polymer results. We also consider the case 



o 

QT 1 of stretching a uniformly charged biopolymer by an electric field, for which we compute the average 



extension and the average shape, which is shown to be trumpet-like. Our results also apply to long 
biopolymers when thermal fluctuations have been smoothed out by a large applied field or force. 
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The mechanical properties of semi-flexible biopolymers are important for their biological 
function. For example, DNA is tightly packed in eukaryotic chromosomes and in viruses, 
while actin and microtubules provide the scaffolding and structure for most animal cells 
In these cases the length scales at which the polymer properties of these macromolecules 
are of biological interest are comparable to their persistence length. Namely DNA packing 
typically involves loops of diameter less than its persistence length of 50 nm, and so does 
looping induced by DNA bound proteins such as the lac repressor, while actin is present in 
cells in the form of bundles and networks in which the typical length of the participating 
polymers is shorter than its 15 /im persistence length. Stretching short strands of DNA is 
also relevant for single molecule experiments involving tethered molecules 

Mechanical properties of semi-flexible polymers are well described by the worm-like chain 

model sl, 4], which treats the polymer as a space curve with a bending energy quadratic 

in the curvature. On the basis of this model, the extension of a molecule in response to 

ft 

an applied force can be calculated [4J, and has been shown to agree with experiments to a 
high level of accuracy. Typically, the method of solution has been to map the calculation 
of the partition function of the worm-like chain to solving a differential equation, either 
numerically Q], or analytically as in the two-dimensional case jsj]. Since experiments have 
usually probed the force response of long DNA molecules, theoretical treatments in the past 
have typically taken the long chain limit, in which case boundary conditions play no role. 
However this is not appropriate for molecules whose length is of the order of a persistence 
length. In fact, the discrepancy between the long chain results and the behavior of short 
molecules is glaringly obvious when we look at the force extension curves themselves. All 
force extension curves for long polymers pass through the origin in the limit as force drops 
to zero, but it has been known for a long time that the short worm-like chain has a finite 
extension at zero force 6|. 

Since the use of short DNA strands as a "force ruler" is becoming routine in force spec- 
troscopy experiments, understanding the entropic elasticity of such strands is of consider- 
able experimental importance . For example, using the formula appropriate for long 
molecules while fitting force-extension measurements done on short molecules leads to incor- 
rect estimates of the persistence length, that could be 2-5 times smaller than the accepted 
value [7|. 

It is interesting to note that similar considerations apply even to long molecules, when 



the molecule in question is being stretched by large forces, assuming that the entropic limit 
is maintained and structural changes to the molecule are not induced. For example, if a 
long molecule is stretched while keeping the tangents at the two ends perpendicular to the 
direction in which the force is exerted, the molecule will bend away from the direction of the 
force at the ends. The stored length in these bends will make the force-extension behavior 
of this molecule different from when it is attached with the tangents parallel to the direction 
of the applied force. It has been proposed that such a mechanism is responsible for the 
observed deviations from ideal force-extension behavior when polymers are stretched by an 



AFM tip 



The short chain and the large force limit discussed above are both characterized by a 
small probability of large deviations from the straight polymer configuration. Throughout 
this paper we refer to this situation as the "fluctuating rod" limit of a semiflexible polymer. 
In this limit we need only take into account small, quadratic fluctuations around the energy- 
minimizing configuration. This leads to a considerable simplification of the mathematics of 
the worm-like chain, and allows for some elegant analytical results, as we demonstrate below. 
Note that an alternative numerical treatment of small-chain entropic properties is provided 
by accounting for the appropriate boundary conditions and finite chain length in the series 
expansion of the partition function for the worm like chain 7L 

Previously, short chains were explicitly discussed in ref. 8|, who also used a harmonic 
approximation to obtain the probability distribution function of the end-to-end vector for 
free polymers in a series expansion. An interesting contribution was made by ref. 9|, who 
discussed the difference in the applicability of thermodynamics to short and long polymers, 
and derived a formula for stretching short polymers. Ref. 10J, along with Monte Carlo 
results, rederived this result for stretching short polymers, and also calculated distribution 
functions and the effect of walls. The effect of stored length due to boundary conditions on 
long polymers stretched by large forces was calculated by ref. jl^], who also considered the 
effect of internal loops on force-extension relations in this regime. In this paper we recover 
prior results for stretching a short polymer by an applied force, using the generating func- 
tional method. We find this field-theoretic method particularly well suited for calculating 
statistical properties of fluctuating rods. To demonstrate this we derive a number of new 
formulas, including the hitherto unsolved problem of the average extension of a charged 
polymer in an electric field. 
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FIG. 1: A cartoon of different ways to stretch single molecules, (a) Magnetic Tweezer apparatus. 
The bead is paramagnetic, and the force exerted is proportional to the spatial gradient of the 
magnetic field, (b) Stretching using a micropipette. The tip is imaged, and its bending is pre- 
calibrated with the corresponding force exerted, (c) Laser tweezers , and (d) stretching using a 
uniform electric field. 



This paper is organized as follows. In SecUwe present the general formalism for calcu- 
lating the force-extension curves and rms-fluctuations of a fluctuating rod under tension, 
using the generating functional method. In SecJTT] we apply this formalism to a polymer 
stretched by a constant force. However single molecule experiments that stretch semiflexible 
polymers by laser tweezers, magnetic tweezers, or a micropipette [4] differ in the way they 
treat the ends of the molecule, which may be free or constrained in several ways (Figfl]). 
These experimental conditions affect the entropy of the molecule, and thereby the force- 
extension relation. We demonstrate this for fluctuating rods by showing that the effect of 
axis-clamping at the two ends of the molecule leads to a measurable effect on force-extension 
curves. 

In Sec JIIII we show that a uniformly charged polymer in a constant electric field behaves 
as if it is being stretched by a force that varies linearly along the contour. Calculating 
the force-extension relationship in this case, using the standard method of obtaining the 
partition function by solving the partial differential equation it satisfies, is not practical, 
since the differential equation in this case is not separable. Previous work in this area has 



relied on either phenomenological arguments 



15], or on linear response approximations 



for weak fields 



161 ] . However, the general method we outline in Sec JT] can easily handle this 



situation as well, and we present new analytical results for the force-extension relation and 
the shape of a uniformly charged fluctuating-rod polymer stretched by an electric field. The 



151 ] . Our results are applicable to 



calculated average shape is in accord with observations 
stiff polymers like actin for a wide range of field strengths, and can be used to measure the 
effective charge density of actin for different salt conditions. 



I. STATISTICAL MECHANICS OF FLUCTUATING RODS 

Since stretching experiments can be performed in both two and three dimensions, we work 
in general ^-dimensions. The objective is to calculate the average extension of the polymer 
in the direction of the applied force or field. We start by making the small-fluctuation 
approximation, which transforms the partition function of the worm-like-chain into a Gaus- 
sian path integral. Using this partition function we obtain the tangent-tangent correlation 
function, from which we derive the average extension of the molecule. 

The Hamiltonian of a worm-like-chain polymer of contour length L, which is stretched 
by a constant force F is, 

™i=f ds {!(l) 2 -™}' (1) 

where (3 = 1 /ksT, s is the arc length and £ is the bending stiffness and is simply related to the 
persistence length, l p — 2£/(d — 1), in (i-dimensions. For double-stranded, DNA £ ~ 50 nm 
and for an actin filament, £ ~ 15 /im. The unit tangent vector t(s) specifies the conformation 
of the chain. The polymer is stretched by tethering the s = end to a fixed support, usually 
a bead held by a micropipette, and pulling on the other end. The usual procedure is to 
attach either a magnetic or a polystyrene bead to the other end of the polymer, and to 
exert a constant force F using magnetic tweezers, optical tweezers, or another micropipette 
(FigHJ). We assume that F is along the e± direction in ci-dimensional space. The components 
of the tangent vector are t(s) = {i-^/l — X^=2^( s ) 2 > ^a( s )> " ' > td{ s )}- In the limit of small 
fluctuations, the tangent vector makes only very small deviations away from the direction 
of the applied force, so that its coordinate along the e.\ direction can be approximated as, 
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In this limit, (dt/ds) 2 ~ "Yli^idti/ds) 2 . Inserting these expressions in the Hamiltonian, 
Eq.flU), and integrating the bending energy term by parts yields, 
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i=2 



dst i (s)0(s)t i (s) + Bi(L, 0) 



(3) 



where we have left out an unimportant constant term. Here the term Bi(L, 0) = — |L 
depends on the boundary conditions, while 



0(s) = -^ + (3F(s) 



(4) 



is a differential operator. We consider the general case of an s-dependent force, F(s), which 
includes the case of stretching by an electric field discussed in Sec JIIIl 

The boundary term Bi(L, 0) depends in general on the tangent vectors at the two ends of 
the polymer, and for the cases discussed here it vanishes. Namely, the tangent vector at the 
end is either unconstrained, in which case dt/ds = for the end-tangent vector, or aligned 
with the force, and then ti = for i = 2,3.... 

The partition function of the fluctuating rod, Z, is a path integral of the Boltzmann factor, 
exp(—/3H[ti(s)]), over all possible conformations of the polymer, where H[ti(s)] is given by 
Eq.Q with the above boundary conditions. To compute the tangent-tangent correlation 
function we employ the generating functional, Z(Ji), which is obtained by adding a source 



term, Jj(s) to the Hamiltonian: 



V[{t(s)}} exp [-/3H[{U(s)}] + Ji{s)ti(s)] . 



(5) 



The correlation function, {U(s)ti(s')) is then obtained by taking functional derivatives 171 ]. 

5 2 log[£(J,)] 



(ti(s)ti(s')) = lim 



G(s,a') 



A^o 5Ji(s)8Ji(s') 

where, G(s, s') is the Green's function of the operator 0(s) defined by, 



(6) 



0(s)G(s,s') 



G(s,s') = 5(s-s') . 



(7) 



Note that, by symmetry considerations, fluctuations transverse to the applied force are 
equivalent for all transverse directions, hence we need to calculate this Green's function for 
only one of the d—1 directions. 



To calculate the average extension in the Xi-direction at any contour position s, we now 
make use of the relation 

Ms)) = Jm(i - {Wf)) , (8) 

which is valid in the small fluctuation approximation regime given by Eq.(j2J). The end-to- 
end extension is X = (xi(L)). The average of each of the other d — 1 orthogonal coordinates, 
Xi (i — 2, . . . , d), is zero, while their root-mean-square (RMS) value is 

(x^s) 2 ) 1 * = (J'ds! J'ds!' (U(s%(s")f) 2 . (9) 

A parametric plot of the mean extension against the root-mean-square value of Xi(s), (i = 
2, . . . , d) using Eq.flS]) and Eq.flH} describes the average shape of the polymer. 

An important question is, under what conditions does the fluctuating rod approximation 
apply? For this we employ the tangent-tangent correlation function which yields a self- 
consistency condition for the small-fluctuation assumption. In two dimensions, the small 
fluctuation approximation, Eq.(j2J), implies that t x « 1 — ty/2 for all s, which is accurate 
to within 1% when ^ < 1/2. In ^-dimensions this suggests the condition that the sum of 
all mean-squared fluctuations at any contour position is at most 1/2. In other words, we 
require that 

(U(s) 2 ) < l/(2d-2) , (10) 

for all s. Below we repeatedly make use of this condition in order to estimate the range of 
experimental parameters (polymer length, magnitude of force, etc.) for which the fluctuating 
rod model is applicable. 

II. STRETCHING BY A CONSTANT FORCE 

To illustrate the method described above, we calculate explicitly the force-extension re- 
lation and the rms-fluctuations of a fluctuating rod stretched by a force applied at one of its 
ends. The experimental setup is as shown in FigfTji. The correlation function, (t(s)t(s')) is 
given by the Green's function of the differential operator, Eq.(j3]). Note that the assumption 
that the tangent vectors at the two ends of the polymer are both aligned with the direction 
of force makes the Green's function vanish at the boundaries. 



The differential equation is solved separately in the domains s < s' and s > s', in which 
the delta function vanishes, and then the solutions are matched at the boundary s = s'. 
This yields, 



G(s, s') 



Af sinh[/«Fs] sinh[/€F(L — s')] = G < (s, s'), when s < s' 
Af' smh[n F (L — s)] smh[n F s'} = G > (s, s'), when s > s' 



where 



(11) 



(12) 



£k B T ' 

is a force-dependent inverse length. Requiring the continuity of G(s, s') at s — s', we find 
Af = Af'. Also, by integrating both sides of Eq.fjTJ) over the interval s 6 (s — e, s + e) and 
taking the limit e — > gives the final boundary condition on G(s, s'): 
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(13) 



Putting it all together, we get, 



sinhf/t^s] smh[n F (L — s') 



k f £ sinhf^L] 
sinh[Kp(L — s)] sinh^s'] 



when s < s' 
when s > s' 



(14) 



G(s,s') 

sinh [KpL] 

This Green's function is plotted in FigO 

When s = s', the correlation function yields the mean-square fluctuations of the tangent 
vector coordinate, (ti(s) 2 ). Plugging this into Eq.©, we obtain 

(xi(s)) s . 2nps cosh [npL] — sinh [upL] + sinh [kf(L — 2s)} . . 

— — = --[fi-\) . (15) 

Now we set s = L in Eq. (fi~5l) to get the force extension relation for a fluctuating rod polymer, 

(16) 



X d — 1 KpL cosh [kfL] — sinh [npL] 

L 2 2k f ^L sinh [KpL] 



A plot of this relation, and its comparison with the long-polymer result is shown in Figj3l 
For d = 3, Eq. lfTB"]) agrees with the force-extension relation computed in Ref. 9[]. 

Experiments that involve stretching short strands of DNA, or other semiflexible molecules, 

fl 

should use Eq.( fl6l) in place of the force-extension relation of Ref. [4|, which is appropriate 
in the long chain limit 181 ] . As we have remarked earlier, use of the long chain formula leads 
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FIG. 2: The Greens function G(s,s'), with up = 1, plotted for s'/L = 0.25 (long dashes), 0.5, 
(solid line) and 0.75 (small dashes). 

to erroneous results when fitting data from stretching short polymers p| . It is reassuring to 
note that, as we shall show below, Eq. (|T6|) reproduces known results in the limits of high 
force and zero force. 

It should be also noted that in the long-polymer limit, the statistical properties of a 
semi-flexible polymer under an applied force depend on one dimensionless ratio, -F£/(fcgT), 
which also delineates the limits of high and low force. In the case of fluctuating rods however, 
both of the two independent dimensionless ratios that can be formed by the three lengths, 
L, £ and present in the functions describing the statistical quantities of interest. 

Also, in this case the high force limit is governed by a different dimensionless number, 
KpL = ^/FL 2 /ik B T. 

We can also calculate the root-mean-square fluctuations of the polymer in the transverse 

directions, using Eq.flU]), 

l I 
{xi(s) ) 2 i 2kfS sinhfK^L] — 3cosh[/ti?L] + 4cosh[/tp(L — s)] — cos!i[kf(L — 2s)] y 

L = V 2K F L 2 £smh[n F L] ) ' 

(17) 

for the transverse directions % = 2, . . . d. A parametric plot of (xi(s)) given by Eq. fll5j) 
against a/ (xi(s) 2 } i = 2, . . . d defines the average shape of the polymer. 

The self-consistency assumption, Eq. fllOj) . is a necessary condition for the validity of 
Eq. (jT6l) and Eq.flTTJ). It is clear from Eq. ffT^j) that the maximum fluctuations of the tangent 




FIG. 3: The force-extension curve for a 100 nm strand of DNA, plotted against force in piconewtons 
(solid line) using Eq. (|16p . The dashed line shows the force-extension curve for a long polymer, 
plotted using the approximate interpolation formula of Ref. 4] 



vector are at the center of the chain, s — L/2. We therefore require that 

2^ 



<^) 2 > < ( M \ 



< 



i 



2(d-l) ' 

Taking d = 2 for illustration, this yields the condition that 



(18) 



2k f £ 



tanh 



KpL 



< 



As the hyperbolic tangent is never greater than one, Eq. (jl"9]) always holds when 



(19) 



(20) 



However, even when the force is weak, Eq.( jl"6l) can apply, provided the polymer is short 



enough compared to its persistence length. In other words, even when 



> 1, the small 



fluctuation assumption can be satisfied, provided the hyperbolic tangent is small enough. 
In this limit Eq. (1191) can be rewritten as 



L 1 , 

7" — 7 l0 § 



(21) 



1 - K F £_ 

Conditions (!20l and (pljl are summarized graphically by the shaded region in FigJU 

It is an interesting exercise to examine the limits of Eq. (1161) in the two cases above, 
described by Eq. (l2"0"l) and Eq. (l2~Tl) . When the polymer is long, so that L/£ ^> 1, the strong 




FIG. 4: The set of parameter values that satisfy the small fluctuation approximation is shown by 
the shaded area. The x-axis is the reduced force, n 2 F ^ 2 = F^/ksT. The y-axis is the polymer 
length in terms of the 3d-persistence length, L/£. 



force condition, Eq.(l20l). ensures that the product kfL ^> 1. Then the force-extension curve 
in Eq. (jl"6l) reduces to 



which, for d = 3, is the well known result for high force stretching [4j. Our result in Eq. (116p 
shows, in fact, that even for a short chain satisfying L/£ < 4, as long as the force is strong 
enough, such that the product KpL ^> 1, we recover a force-extension relation of the form 
Eq.©. 

The other limit is the case when the force is weak such that KpL < 1, then the force 
extension curve in Eq. (jT6l) becomes 

T^-^fKi)-^© 3 ) • (23) 

This is different from the behavior of long worm-like chain polymers under weak forces, 
in which case the force-extension relation is proportional to the force |4j, with no constant 
term. Note also that when the force goes to zero, the relative extension of the polymer tends 
to 1 — which is a result one finds in Landau and Lifshitz for d = 3 [6]. 

What about the case when 1/(2k_f£) > 1? Then we must have tanh(K^L/2) < k_f£, or 



in other words, to a first approximation, 



L < 2£ , (24) 

which can be confirmed by a look at FigOJ Hence polymers up to one persistence length in 
two dimensions, or two persistence lengths in three dimensions, can be regarded as satisfying 
the small fluctuation approximation for any applied force. 



A. Effects of Axis Clamping 

We now show how axis clamping affects the force extension for a fluctuating rod. If the 
free end of the polymer is held rigidly, such that the end-to-end separation vector of the 
polymer is constrained to be collinear with the direction of the force, there is an additional 
constraint, 

dsU(s) = Q for i = 2,3,...,d. (25) 

Note that this constraint depends upon the entire conformation of the chain, hence it is not 
a priori clear when it can be ignored. 

We can compute the correlation function, (tj(s)ij(s')), for all i = 2, 3, . . . , d, by forcing 
the new constraint in Eq.( !25l) using the Dirac delta function, 



(ti(s)ti(s')) 



Jvu{s) ti{s)U(J) s(J dsU{s) 


exp 


Jds^-^tiWOWUis^+pFL 


Jvti(s)5^J dsti(s)jexp 


£ds ^ — t^OisMs)^ + (3FL 





(26) 



We use the integral representation of the delta function in Fourier space, add the source 
terms as before and compute the Gaussian integrals by completing the square. This gives 
us an additional term in the correlation function, 



(t i (s)t i (s')) = G{s,s') 



J dsiG(s,si^j(J ds 2 G(s', s 2 ] 



(27) 



ds x \ ds 2 G(s 1 ,s 2 ) 
Jo 

where G(s,s') is the Green's function of the operator 0(s) shown in Eq. (114ft . Now using 



Eq. (j2"71) in Eq.(jSD, we obtain the average extension, 



X (d — l)(npL cosh[/tpL] — sinhf/tpL]) 

L 4K F £Lsmh[K F L] 
{d — 1)(kfL cosYl[kfL\ — 3sinh[KpL] + 2npL) 

(28) 



2n F £L 2 + 2np£L 2 coshf/tpL] - 4k 2 f £L sinh[K F L] 

The first two terms in Eq. fl28|) are identical to Eq. f|T6l) derived above. Therefore, we find 
that the loss in entropy due to axis clamping of the two ends of the polymer, leads to a 
small but measurable correction to the average extension. This is plotted in Figj5]for d = 2. 
We can also use Eq.flHJ) and calculate the average shape of the polymer, defined by Eq.flSJ) 



and Eq.flUJ), by plotting (xi(s)) against a/ (x2(s) 2 } parametrically, as shown in Figl5b. 

We calculate the range of validity of Eq. (l28]) by imposing Eq. ffTO]) as before. The param- 
eter values that satisfy this condition are displayed in the inset of Figj5l It can be shown 
that when L/£ < 6, Eq. fl28l holds for all F, while for long polymers axis clamping has no 
effect on the force-extension relationship, and for kfL > 1 we again recover the well-known 
result 4], Eq.(F22]). 

We can expand Eq. (1281) in the limit of small forces to obtain an equation analogous to 
Eq. (1231) for the case of axis clamping. Interestingly, for d = 3 we obtain, 

* = 1 -±(k) + (L) 3 . (29) 

L 30 \CJ 25200k B T\£J K ' 

Note that as F — > 0, X/L approaches 1 — L/(30^), while Eq.( l23l shows that without the 
entropic constraint, X/L approaches 1 — L/(6£) in the limit of zero force. Hence we see a 
different behavior at low forces due to the entropic constraint at the end. 

An experimentally relevant situation is illustrated in Fig. [Tb, when the polymer is 
stretched using a micropipette and the force exerted is calculated by the bending of the 
micropipette tip. An equivalent setup uses a bead attached to a glass fiber and measures 
the bending of the fiber. These methods have been used earlier, for example to stretch 
11 ] and long DNA molecules Ijj]. Stretching using an AFM also involves axis 



chromosomes 



clamping, though here the applicability of our formulae depends on the other boundary con- 
dition, i.e. the collinearity of the first and last tangent vectors with the force, being satisfied. 
In principle, the entropic effect of axis-clamping may also be observed in stretching with a 
laser tweezer since the optical bead is held in a three dimensional trap formed by the laser 
beam (FigJTb). In practice however the trap is not a perfect clamp. First, it is significantly 




FIG. 5: (a) The relative extension of a polymer X/L pulled by a constant force in the presence of 
axis-clamping is plotted against the reduced force F^/(ksT), for a short polymer (solid line) for 
d = 2. The dot-dashed line is the extension without the axis-clamping condition, Eq. (|16p while 
the dashed curve is the long chain result [5|. The reduced polymer length is L/£ = 0.4 (~ 20 nm 
for DNA and ~ 6/^m for actin). In terms of force, F£/(kBT) = 1 means 0.082 pN for DNA, and 
0.27 fN for actin. Inset shows the range of allowed parameter values (shaded area), (b) The shape 
of a polymer stretched with a laser tweezer is plotted parametrically in units of polymer contour 
length L for a short polymer (L/£ = 0.4). The solid curve is for F£/(kBT) = 1, and the dotted 
curve is for F£/(k B T) = 1000. 

weaker in the z-direction, so some fluctuations do occur. Second, it is very hard to control 
where on the surface of the bead the polymer binds, so collinearity may be only approxi- 
mately satisfied, since the bead is about three orders of magnitude larger that the polymer 
thickness. Finally, the bead is free to rotate in the trap, to the extent allowed by the twist 
elasticity of the polymer 131 ] . Experimentally the effect of axis clamping might be relevant 



for stretching very short DNA strands - a O.lpN force will stretch a 300nm strand of DNA 
to 0.86 of its contour length if stretched with axis-clamping, compared to 0.81 without it. 
The difference might also be observable when pulling on actin of about a persistence length 
in size, but in the range of femtonewton forces. 



III. FLUCTUATING ROD IN AN ELECTRIC FIELD 



Next we consider stretching of a charged polymer by a constant electric field. Here, one 
end of the polymer is tethered, and when the field is switched on, the polymer is observed 



to extend 



15 



20| | . Because the polymer is in a charged aqueous environment, the molecular 



mechanism of the electrophoretic stretch is complicated 2Qj,|22j. Here we abstract from these 
difficulties and assume that the polymer responds to the field £ as a uniformly charged rod 
with charge per unit length A. 

As in the previous section the contour length along the polymer is denoted by s, with 
s = denoting the tethered end. If the position vector of the segment s along the contour 
is denoted by r(s), the interaction potential of the polymer with the field can be written as, 

Hi[r(s)] = -XE- I ds (r(s) - r(0)) . (30) 
Jo 

Note that we are not only assuming that the effective charge density is constant, but that 
it also remains unchanged during the course of stretching. We can express Eq. (1301) in terms 
of tangent vectors, 

H I [t(s)] = -\E- [ ds I ds't(s') (31) 



JO JO 

We now change the order of integration. Instead of integrating first over s' from to s, 
and then over s from to L, we can equivalently first integrate over s, from s' to L, then 
integrate over s', from to L. Hence we get, 

Hi = -XE ■ ds dst(s') 

J0 J s' 

= - [ ds'XE(L - s') ■ t(s') , (32) 
Jo 

Thus a constant electric field stretches on a uniformly charged polymer as though it were 
subjected to a contour dependent force, F(s) = EX(L — s). This force is zero at the free 
end, s = L, and reaches a maximum of \E\XL at the tethered end. 



The s-dependent potential Hj makes the application of spectral methods to evaluating the 
partition function not practical. Namely, in this case the calculation of the eigenf unctions 
for the spectral representation of the partition function maps to a Schrodinger-like equation 
with a time dependent potential, making the analysis quite complicated. The generating 
functional formalism, on the other hand, can be employed without much difficulty 

The boundary conditions appropriate for this situation are £j(0) = 0, i — 2,3, ... ,d at the 
tethered end, and ^4^-\ s =l = 0, for all i at the free end. The boundary term in Eq.flS]) once 
again vanishes. The correlation function (9i(s)9i(s')) is now the Greens function, G(s,s'), 
of the operator 

0{s) = + PXE(L - s) , (33) 

with the boundary conditions above. 

The solution to Eq.Q that satisfies the boundary conditions is calculated in the same 
way as for the constant force case discussed above, but now in terms of Airy functions, Ai 
and Bi, and their derivatives, Ai' and Bi', 



, t) = 7T (Ai[ Kg (L - s)]Bi'[0] - Ai'[0]BiML - s)]) (Ai[K E (L - s')]Bi[ KE L] - Ai[ KE L]Bi[ KE (L - s')}) 

(34) 

Eq. (134D is the solution for s > s', while the case s < s' follows from symmetry. 
The relevant, electric-field dependent length scale is given by k e 1 , where 

K E = (PXE/S)*. (35) 

It is interesting to note that this length scales as i? -1 / 3 , which is different than the —1/2 
power that characterizes the dependence of k f 1 on the force. This could not have been 
predicted by dimensional arguments. In fact, since EXL is a force, one might have guessed 
the same scaling for fOr Kn . 

We can now calculate the relative extension of the polymer in the direction of the field 
using Eq.©: 

X = (d - 1) (Ai'[K g L]Bi'[0] - Ai'[0]Bi'[/^L]) 

L 2K 2 E ^L(Ai[K E L]Bi'[0]- Ai'[0]Bi[K E L]) ' 1 ' 

Eq. (|36l) is one of the main new results of this paper. It provides, for the first time, an exact 

expression for the field-extension relation for a charged semi-flexible polymer in an electric 

field. 



Since the quadratic Hamiltonian, Eq.(j3]), implicitly assumes that the fluctuations of the 
tangent vector are small, we require (ti(s) 2 ) to be small for our analysis to be self consistent. 
Once again we impose the self-consistency condition, Eq. fllOp . on the variance of the fluctu- 
ations, (ti(s) 2 ). It is clear from the boundary conditions, that the tangent vector is fixed at 
the s = end and completely free at the s = L end, that the maximum fluctuations of the 
tangent vector occurs at s — L. Thus, the self-consistency requirement can be written as 

= Ai[0]Bi[^]-Ai[^L]Bi[0] < 1 
x ' 0i (Ai[K B L]Bi'[0] - Ai'[0]Bi[K E L]) " 2(d — 1) 

When the argument of Airy function is large, Ai converges to whereas Bi diverges to +oo 
exponentially as 



2, 3 n 
'3' 



Ai[z] 



m[z ] ~ (3g) 



exp[g^2 



,'4 



Thus, when k e L ^> 1, the condition in Eq. (l37|) reduces to 

WW = -A £ WT) (39) 

Again we get a ratio of length scales, i.e. the strong field condition is equivalent to 

k e £ > 2.7A(d - 1) , (40) 

where we have used Ai[0]/Ai'[0] ~ 1.37. The small fluctuation approximation therefore 
requires that the length-scale, n E x be less than about a third of the persistence length. 

On the other hand, when k e L <C 1, we can Taylor expand Ai[keL\ around as Ai[K^L] ~ 
Ai[0] + K^LAi'fO] and similarly for Bi. Then, Eq. fl37|) becomes 

J^W^T)- (41) 

Thus Eq. fl36|) is valid for all field strengths for polymers less than £/2 in two dimensions, and 
£/4 in three dimensions. For longer polymers, stronger fields are needed for its applicability. 
For example, Eq. fl36|) is applicable to a molecule of actin that is about 15/im in size, and 
is being stretched by an electric field of at least 0.02V/cm in strength. But it is not useful 
for a molecule of DNA longer than about 25nm; the minimum field required for a 50nm 
DNA molecule for the small fluctuation approximation to be appropriate is of the order of 



10 2 V/cm. The inset of Fig. O shows the region of parameter space, where in d = 2 Eq. (j3"7|) 
holds for the whole range of KgL. 

When the chain is long and the applied field is strong such that E ^> (/5A^ 2 ) -1 , the 
expression of extension reduces to 



It should be noted that it is possible to derive this equation using approximate phenomeno- 

nn 

logical arguments [4J, [15( . A comparison with the constant force case here is very interesting. 
As we have shown, the length scale, k f 1 in the constant force case scales as F~ lj/2 in contrast 
with k^ 1 that scales as E~ 1 ^. In the high field limit however, the dependence of the relative 
extension on the electric field becomes E" 1 ^ 2 , just as the dependence of extension on force 
in the high force limit. 

In the other limit, when the polymer is short so that L/£ < l/(2d — 2), and the electric 
field is small, E <C (/3A£ 2 ) -1 , Eq. (l36|) yields the limiting linear response, 

L 4 V?/ 40k B T \U ^ 1 

Eq. (j43D has been derived for d = 3 using different methods in Ref. 16]. 

Again we can compute the average shape, defined as before by Eq.® and Eq.Q, of 
the polymer under stretch. In experiments on long DNA molecules, it has been reported 
that the polymer assumes a trumpet-like shape, reflecting the fact that force is stronger 
towards the grafted end and tends to vanish at the free end jig . The integral in Eq.fjH]) can 
be calculated analytically, while Eq.Q needs to be calculated numerically. A trumpet-like 
shape emerges naturally as the result of our computation, and is shown in Fig. [6j 

In experiments that stretch single molecules by electric fields, the polymer molecule, is 
in an ionic solution, and an electric field will drive a net current The current creates a 
flow in the system, and the net deformation of the molecule, it has been argued, is therefore 



a combination of the flow and the field 



20 



221 ] . In experiments performed with long 



DNA molecules of different lengths, it has been shown that Eq. (|42|) fits the data only with 
different effective charge densities for different lengths of the molecule, thus questioning 
whether the effective charge density is a physically meaningful parameter 20J]. However we 
feel that interpretation of the results have been hampered by the lack of precise formulas 
for the electric field induced stretching case. While this continues to be the case for long 




FIG. 6: (a) Relative extension, X/L is plotted against the reduced electric field, j3XE^ 2 for a 
polymer of size L/£ = 0.4 for d = 2. Inset shows parameter values allowed by self-consistency 
(shaded area). For actin, if A ~ le/nm, fiXEt^ 1 = 1 is about 10 -3 V/cm, for DNA with 0.6e/nm it 
is about 170V/cm. Note that we are assuming that the screened charge density is about one-tenth 
of the bare charge density. This is consistent with experiments on DNA [15]. (b) The shape of a 
short polymer stretched by a constant electric field in two dimensions is plotted parametrically in 
units of polymer contour length L/£ = 0.4. The solid curve is for j3\E£? = 1, and the dotted curve 
is for (3\E£ 2 = 1000. 

DNA molecules, Eq. (1361) provides a means of testing these questions using actin molecules. 
Quantitative comparison of our results with experiments would be relatively easy for actin 
polymers, since the required electric fields for the ffuctuating-rod assumptions to hold are 
small, and easily produced in the laboratory. Assuming that charge screening reduces bare 
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FIG. 7: The relative extension of a 15 micron actin molecule in a constant electric field compared 
to the case when it is stretched by a constant force due to a single charge on the free tip of the 
polymer. When this charge is chosen to be exactly 1/4 of the total charge on the polymer, the two 
curves agree well over most of the range of applicability. 

actin charge to about 1 electrons/nm, our results are valid for all actin lengths at electric 
field strengths starting from 0.03 V/cm upwards. From force-extension measurements in such 
experiments, the effective charge density of actin of different lengths and under different ionic 
conditions can be extracted. 

In view of the non-uniform nature of the effective force a charged polymer experiences in 
a constant electric field, it may be useful to ask what is the effective constant force, which 
would lead to the same relative extension as a given electric field strength. Eq. (142!) appears 
to suggest that at high forces the polymer behaves as if it is being stretched by a constant 
force equivalent to that produced by one-fourth of the total charge concentrated at the end 
of the chain. It turns out that over most of the regime of applicability of our result, that is 
a reasonable approximation, as demonstrated in FigJTl 

IV. CONCLUSIONS 

Biopolymers like DNA are fast emerging as tools - rulers and templates - for a number 
of biophysical applications, such as measuring DNA-protein interactions. Some of these uses 
are due to the remarkable fact that a simple physical model, the worm-like chain, explains 



DNA entropic elasticity with great accuracy, and small deviations from this behavior, say 
due to bound proteins, can be detected. However entropic effects are subtle and may be 
affected by experimental conditions in ways that are not obvious. For example it was recently 
pointed out that tethered bead experiments, that use a short DNA strand attached to a bead, 
need to take into account entropic exclusion forces of the bead with the wall [2l| . Similarly, 
experiments involving entropic elasticity need to account for the influence of experimental 
conditions on the entropy of the molecule, that may quantitatively change the force-extension 
relation. 

In this paper we have theoretically examined entropic effects of boundary conditions on 
force-extension curves. We demonstrated that the force-extension response of fluctuating 
rods is qualitatively different from that of long polymers, for which the contour length 
is much greater than its persistence length. We also derived analytic expressions for the 
force-extension relation and root-mean-square fluctuations of such a polymer stretched by 
a constant force. Finally, we showed that in the fluctuating rod limit the constraint of 
axis- clamping, which might be imposed by single-molecule stretching techniques such as 
laser tweezers, can lead to measurable effect on force-extension curves. More importantly, 
blindly employing the long-polymer formulas to data obtained for short chains might lead 
to erroneous conclusions [7]. 

Single molecule experiments using electric fields to stretch molecules may also become 
a useful tool, but progress is hampered by the lack of analytic expressions for the force- 
extension curves in this case. We examine the case of a polymer stretched by an electric 
field, and show that it behaves as if it is being stretched by a non-uniform force, which 
increases linearly from the free end of the chain. We derive an analytic expression for the 
field-extension relation and average shape of the polymer stretched by an electric field, in 
the fluctuating rod limit. 

In all cases considered here, the force extension formulae, Eq. (1281) and Eq. (1361) . are func- 
tions of three length scales, two of which are the contour length L and the 3<i-persistence 



length £. For stretching by a constant force, the third length scale is kJ, 1 = ^ksT^/F, 
while for a constant electric field it is = {ksT^/ E\)z . It is noteworthy that these 
length scales do not follow from dimensional arguments alone due to the presence of di- 
mensionless combinations of polymer parameters. Since these length scales are related to 
the decay of fluctuations in the directions transverse to the force, and become smaller as 



the polymer becomes more straight, they may be interpreted as confinement length-scales, 
pointing towards deeper analogies between stretched and confined polymers |23J that warrant 
further exploration. 
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